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Q ■ Abstract. 

T _j_ , We develop a theory for option pricing with perfect hedging in an inefficient market 

model where the underlying price variations are autocorrelated over a time r > 0. This 
■ is accomplished by assuming that the underlying noise in the system is derived by an 

Ornstein-Uhlenbeck, rather than from a Wiener process. After obtaining an effective 



onc-dimcnsional market model, we achieve a closed expression for the European call 
■ price within the Black-Scholes framework and find that our price is always lower than 

the Black-Scholes price. We obtain the same price if we start from a modified portfolio 
• although now we get a different hedging strategy than that of Black-Scholes. We 

compare these strategies and study the sensitivity of the call price to several parameters 
where the correlation effects are also observed. 

s 

C , 1. Introduction 

o 
o 

Fischer Black and Myron Scholes (1973) and Robert Merton (1973a) obtained a 
fair option price assuming severe and strict theoretical conditions for the market 
behavior. The requirements under which these were developed include: (i) Absence 
of arbitrage opportunities, i.e., identical cashflows have identical values (Sharpe (1964); 
Cox and Ross (1976)). (ii) Efficient market hypothesis, i.e., the market incorporates 
instantaneously any information concerning future market evolution (Fama (1965)). (iii) 
Existence of a unique riskless strategy for a portfolio in a complete market (Markowitz 
(1952)). Due to the random character of stock market prices (Cootner (1964)), the 
implementation of these conditions, especially condition (ii), indicates that speculative 
prices are driven by white (i.e., delta-correlated) random processes. At this point, one 
has to choose between a Gaussian white process (Black and Scholes (1973); Merton 
(1973a)) or a white jump process. In this latter case and due to requirement (iii), the 
jump lengths also have to be known and fixed (Merton (1976)). There are no other 
choices for modelling market evolution if the above requirements and ideal conditions 
are to be obeyed (Cox and Ross (1976)). 

From these three assumptions, condition (ii) is perhaps the most restrictive and, 
in fact, disagrees with empirical evidence since real markets are not efficient, at least at 
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short times (Grossman and Stiglitz (1980); Fama (1991)). Indeed, market efficiency is 
closely related to the assumption of totally uncorrected price variations (white noise). 
But white noise is only an idealization since, in practice, no actual random process 
is completely white. For this reason, white processes are convenient mathematical 
objects valid only when the observation time is much larger than the autocorrelation 
time of the process||. And, analogously, the efficient market hypothesis is again a 
convenient assumption when the observation time is much larger than time spans in 
which "inefficiencies" (i.e., correlations, delays, etc..) occur. 

Alternative models for describing empirical results of the market evolution have 
been suggested (Mandelbrot (1963); Fama (1963)). In each of these, an option price 
can be obtained only by relaxing some or even all of the initial Black-Scholes (B-S) 
assumptions (Figlewski (1989); Aurell et al (2000)). Our main purpose in this paper 
is to derive a nontrivial option price by relaxing the efficient market hypothesis and 
allowing for a finite, non-zero, correlation time of the underlying noise process. As a 
model for the evolution of the market we choose the Ornstein-Uhlenbeck (O-U) process 
(Uhlenbeck and Ornstein (1930)) for three reasons: (a) O-U noise is still a Gaussian 
random process with an arbitrary correlation time r and it has the property that when 
t = the process becomes Gaussian white noise, as in the original Black-Scholes option 
case, (b) The O-U process is, by virtue of Doob's theorem, the only Gaussian random 
process which is simultaneously Markovian and stationary (Doob (1942)). In this sense 
the O-U process is the simplest generalization of Gaussian white-noise, (c) As we will 
see later on, the variance of random processes driven by O-U noise seems to agree with 
the evolution of market variance, at least in some particular but relevant cases. 

The Ornstein-Uhlenbeck process is not a newcomer in mathematical finance. For 
instance, it has already been proposed as a model for stochastic volatility!] (SV). Our 
case here is rather different since, contrary to SV models, we only have one source of 
noise. We therefore suggest the O-U process as the driving noise for the underlying price 
dynamics when the volatility is still a deterministic quantity (Dumas et al (1998)). 

The autocorrelation in the underlying driving noise is closely related to the 
predictability of asset returns, of which there seems to be ample evidence (Breen and 
Jagannathan (1989); Campbell and Hamao (1992)). Indeed, if for some particular stock 
the price variations are correlated during some time r, then the price at time ti will be 
related to the price at a previous time t\ as long as the time span t 2 — t x is not too long 
compared to the correlation time r. Hence correlation implies partial predictability. 
Other approaches to option pricing with predictable asset returns are based under the 
assumption the market is still driven by white noise and predictability is induced by the 
drift (Lo and Wang (1995)). Since the B-S formula is independent of the drift, these 
approaches apply B-S theory with a conveniently modified volatility. Our approach here 
is rather different because we assume the asset price variations driven by correlated noise 

§ Throughout this paper we will use the terms "correlation" and "autocorrelation" without distinction. 
|| Hull and White (1987); Scott (1987); Stein and Stein (1991); Heston (1993); Ghysels et al (1996); 
Heston and Nandi (2000). 
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-which implies some degree of predictability. 

Summarizing, our purpose is to study option pricing and hedging in a more realistic 
framework that of white noise process presented by Black and Scholes. Our model 
includes colored noise and the dependence of the volatility on time. Both are empirically 
observed in real markets (Bouchaud and Potters (2000)). Empirical characteristic time 
scales are at least of the order of minutes and can affect option prices particularly when 
the exercising date is near and speculative fluctuations are more important. Presumably, 
this effect is negligible when correlation times are shorter (much shorter than time to 
expiration). In any case, it is interesting to know how, and by how much, the option price 
and its properties are modified when correlations in the underlying noise are significant. 

The shortest way of getting the call price, and hence quantifying the effect of 
correlations on prices is by martingale methods. Unfortunately, this procedure does not 
guarantee that we obtain the fairest price since arbitrage and hedging are not included 
in this approach. It is therefore our main objective to generalize B-S theory not only 
to get a new call price but, more importantly, to obtain a hedging strategy that avoids 
risk and arbitrage opportunities. 

From a technical point of view, we apply the B-S option pricing method after 
projecting the two-dimensional O-U process onto a one- dimensional diffusion process 
with time varying volatility. As we will show, this projection allows us to maintain 
the conditions of a perfect hedging and the absence of arbitrage. Moreover, the price 
obtained using this way completely agrees with the price obtained using two alternative 
and different methods. One of them is based on martingale theory, and the other one 
develops a new option pricing with a modified portfolio containing secondary options 
instead of the underlying stock. 

The paper is divided into eight sections. In Section |2| we present our two- 
dimensional stochastic model for the underlying asset. In Section |3] we find the O-U 
projection onto the stock price correlated process. Section £| concentrates on the B- 
S option price derivation with the projected process, and Sections [| and |^ show the 
consistency of this derivation by using two alternative methods for obtaining the option 
price. The greeks and the new hedging are presented in Section [7]. Conclusions are 
drawn in Section |8] and technical details are left to the appendices. 



2. The asset model 



The standard assumption in option pricing theory is to assume that the underlying price 
S(t) can be modelled as a one-dimensional diffusion process: 

^ = f,dt + adW(t), (1) 

where W(t) is the Wiener process. In the original B-S theory both drift [i and volatility 
a are constants. Other models take /i = /i(t, S) and a = a(t, S) as functions of time 
and underlying price (Cox and Ross (1976); Bergman et al (1987)). The parameter a 
is assumed to be a random quantity in the SV models. 
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Notice that if the time evolution of the underlying price is governed by Eq. (Jl|) then 
S(t) is an uncorrelated random process in the sense that its zero-mean return rate defined 
by Z(t) = d\n S / dt — ft is driven by white noise, i.e., E[Z{ti)Z{t2)} = cr 2 5(ti —£2) where 
8(t) is the Dirac delta function^]. Hence, the asset model immediately incorporates price 
return effects and meets the efficient market hypothesis. 

As a first step, we assume that the underlying price is not driven by the Wiener 
process W(t) but by O-U noise V(t). In other words, we say that S(t) obeys a singular 
two-dimensional diffusion 
dS(t) 



S(t) 



fidt + V(t)dt (2) 



dV(t) = dt+- dW{t), (3) 

T T 

where r > is the correlation time. More precisely, V(t) is O-U noise in the stationary 
regime, which is a Gaussian colored noise with zero mean and correlation function: 

E[V( tl )V(t 2 )} = g e -I^IA. (4) 

We call the process defined by Eqs. (0)-([|) singular diffusion because, contrary to 
SV models, the Wiener driving noise W{t) only appears in one of the equations, and 
this results in a singular diffusion matrix (Gardiner (f985)). Observe that we now deal 
with autocorrelated stock prices since the zero-mean return rate Z(t) is colored noise, 
i.e., E\Z (t{) Z (t-i)] = [a 1 1 It) exp[— \t\ — t%\lr\. Note that when r = this correlation 
goes to a 2 5(ti — £2) and we thus recover the one- dimensional diffusion discussed above. 
Therefore the case of positive r is a measure of the inefficiencies of the market. 

There is an alternative, and sometimes more convenient, way of writing the above 
equations using the asset return R(t) defined by 

R(t)=\n[S(t)/S ], 

where So — S(to) and to is the time at which we start observing the process (|2])-(|3]). 
Without loss of generality this time can be set equal to zero (see Appendix A). Instead 
of Eqs. (0)-([3|), we may have 
dR{t) 



dt 



H + V(t) (5) 



dJ ^ = \ [-v® + o-m , (6) 

where £(£) = dW(t)/dt is Gaussian white noise defined as the derivative of the Wiener 
process. This process exists in the sense of generalized random functions (Lighthill 
(1958); Stratonovich (1963)). The combination of relations in Eqs. @ and (^) leads to 
a second-order stochastic differential equation for R(t) 
d 2 R(t) dR(t) 



T- 



dt 2 dt 



% We recall that 6(x) is a generalized function with the properties: S(x) = for 1 / and 
J™ 5(x)dx = 1 (Lighthill (1958); Stratonovich (1963)). 



Option pricing and perfect hedging on correlated stocks 5 

From this equation, we clearly see that when r = we recover the one- dimensional 
diffusion case (|1])0 . We also observe that the O-U process V(t) is the random part of 
the return velocity, dR/dt, and we will often refer to V(t) as the "velocity" of the return 
process R(t). 

In Appendix A, we give explicit expressions for V(t) and for the return R(t). We 
prove there that R(t) is a non- stationary process with the conditional mean value 

m(t,Vo)=E[R(t)\Vo] = f£t + T{l-e-^ T )V Q , (8) 

where V(0) = Vq is the initial velocity. The conditional return variance, 

K n {t)=E[(R(t)-m{t,V ))*\V ], 



is given by 



K n (t) = a 



t - 2r (l - e-*/ T ) + - (l - e 



T f ■ -2t/r 



(9) 



2 

We also give in Appendix A explicit expressions for the joint probability density function 
(pdf) p(R, V, t), the marginal pdf 's p(R, t) and p(V, t) of the second-order process R(t), 
and the marginal pdf p(S,t\So,t ) of the underlying price S(t). We also show that the 
velocity V(t) is, in the stationary regime, distributed according to the normal density: 

Pst(V) = ^L=e-^ 2 . (10) 
Jtco- 2 /t 

Suppose now that the initial velocitiy Vq is random with mean value E[Vq] and 
variance Var[Vb]. Thus, the return unconditional mean and variance read 

E[R(t)} =fjLt + r(l-e- t ^)E\y ], 

E[(R(t) - E(R(t)) 2 ] = K u (t) + r (l - e"*/ T ) Var[F ]. 

If, in addition, we assume that the initial velocity Vq is in the stationary regime then 
E[Vq] = and Var[Vo] = a 2 /2r. In this case, the return unconditional mean value is 

m(t) = E[R(t)] = ^t, 
and the return unconditional variance 

n(t) = E[(R(t) - m(t)) 2 } 
reads (c/. Eq. ©) 

K (t) =a 2 [t- r (l-e-*/ T )] . (11) 
A consequence of Eq. (|TT|) is that, when t <r, the variance behaves as 

K(t) ~ (a 2 /2r)t 2 , (t < r). (12) 
Equation ([TT|) also shows a crossover to ordinary diffusion (B-S case) when t^> t: 

n(t) ~ a 2 t, (t > r). (13) 



+ In the opposite case when r = oo, Eq. (^) shows that = 0. Thus V(t) is a constant, which we 

may equal to zero, and from Eq. (^]) we have S(t) = Soe^'. Therefore, the underlying price evolves as 
a riskless security. Later on we will recover this deterministic case (see, for instance, Eq. (|35|)). 
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Figure 1. The variance of the underlying asset as a function of time (in logarithmic 
scale). Circles correspond to empirical variance of S&P 500 cash index from 1988 to 
1996. Solid line represents the theoretical variance, Eq. (7) with r = 2 minutes. The 
dashed line is the B-S variance a 2 t. In both cases a — 3.69 x lO^min -1 / 2 which 
approximately corresponds to an annual volatility a = 11%. 

In Fig. |I|, we plot n(t) along with the empirical variance from data of the S&P 500 cash 
index during the period January 1988-December 19960. The dashed line represents 
results obtained by assuming normal- diffusion n(t) oc t. Observe that the empirical 
variance is very well fitted by our theoretical variance n(t) for a correlation time r = 2 
minutes. Furthermore, the result of this correlation affects the empirical volatility for 
around 100 minutes. These times are probably too small to affect call price to any 
quantifiable extent. However, the S&P 500 is one of the most liquid, and therefore most 
efficient, markets. Consequently, the effect of correlations in any other less efficient 
market might significantly influence option prices and hedging strategies, and this is the 
main motivation for this work. 

3. The projected process 

One may argue that the O-U process (|||)-(|3|) is an inadequate asset model since the share 
price S(t) given by Eq. (fj) is a continuous random process with bounded variations. 
As Harrison et al (1984) showed, continuous processes with bounded variations allow 

* Tick by tick data on S&P 500 cash index has been provided by The Futures Industry Institute 
(Washington, DC). 
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arbitrage opportunities and this is an undesirable feature for obtaining a fair price. 
Thus, for instance, arbitrage would be possible within a portfolio containing bonds and 
stock whose strategy at time t is buying (or selling) stock shares when fi + V(t) is greater 
(or lower) than the risk- free bond rate (Harrison et al (1984)). 

In our case, however, the problem is that in the practice the return velocity V(t) is 
nontradable and its evolution is ignored. In other words, in real markets the observed 
asset dynamics does not show any trace of the velocity variable§]. This feature allows 
us to perform a projection of the two-dimensional diffusion process [S(t),V(t)} onto a 
one-dimensional equivalent process S(t) independent of the velocity V. We will show in 
this section that the projected process S(t), which is equal to the actual price S(t) in 
mean square sense, obeys the following one-dimensional SDE 

= [// + k(T - t)/2]dt + y/k(T-t)dW{t), (14) 

where n(t) is given by Eq. (|TI|), and the dot denotes time derivative. Therefore, the price 
given by Eq. ( TT4] ) is driven by a noise of unbounded variation, the Wiener process, and 
the Harrison et al (1984) results do not apply. In consequence, the O-U projected process 
is still a suitable starting point for option pricing since it does not permit arbitrage. 



3.1. Derivation of the one- dimensional SDE 



Note that the dynamics of the return R{t) = hi[S(t)/ S ] is given by the second-order 
SDE (0) which includes the stochastic evolution of the velocity V{t). Let us now obtain 
a first-order SDE describing the price dynamics when velocity V(t) has been eliminated. 

The starting point of our derivation is the marginal conditional density 
p(R, t\R , to; Vq). This density is given by Eq. ( |A.14j ) of Appendix A and when t ^ it 



reads 



p(R,t\Ro,t ; V 



2irK 11 (t-t 



cxp 



[R-Rq- m{t - t 



Vo)} 2 



0, 



2K n (t-t 



(15) 



where m(t, Vq) and Ku(t) are given by Eqs. @ and @. Note that 

p(R, t\Ro, to', Vq) is the solution of the following partial differential equation 



dp 



[i + V e~ 



(t-to)/r~ 


dp O 2 


1 _ e -(t-to)/r 




dRo 2 





2 d 2 p 



(16) 



with the final condition p(R,t\R ,t] Vb) = 5(R — Ro). Observe that the Eq. fll6|) is a 
backward Fokker-Planck equation whose drift, [i + Voexp[— (t — t )/T], and diffusion 
coefficient, \a 2 [1 — exp[— (t — to)/r]] 2 , are both functions of t — to. As is well-known, 



Indeed, knowing V(t) would imply knowing the value of the return R(t) at two different times, since 

R(t) - R(t - e) 



V(t) = lim 

£-.0 + 



//.. 



Obviously, this operation is not performed by traders who only manage portfolios at time t based on 
prices at t and not at any earlier time. 
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there exists a direct relation between the Fokker-Planck equation and the SDE governing 
the process (Gardiner (1985)). In our case, the corresponding SDE is 

dR(t Q ) = [fi + V e- {t - tt,)/T ] dt Q + a[l- e ~ {t ~ to)/r ] dW(t ), (17) 

and its formal solution is 

R(t) = R(t ) + n{t - t ) + V r 



1 -e 

* r 



to 



(t-t(,)/r 



3.2. Equality of processes in mean square sense 

To avoid confusion, let R(t) be the solution of the first-order SDE (|T7|), i.e., R(t) is the 
projected process given by Eq. (|18|). And let R(t) be the solution of the second-order 
SDE (0) where the dynamics of the velocity is still taken into account. Thus, R(t) is 
explicitly given by Eq. (|A.1|) of Appendix A. 

We will now prove that R(t) and R(t) are equal in mean square sense. That is: 

E [(R(t) - R(t)) 2 ] = 0, for any time t. (19) 

In effect, from Eq. (|l"8"D and assuming, without loss of generality, that t = and 
R(t ) = we have 

rt 



R(t) =fit + V t(1 - e 



-t/r\ 



1 - e 



-(t-tl)/T 



Z(h)dt 



(20) 



where £(t\) = dW{t\)/dti is the Gaussian white noise. On the other hand, from 
Eq. ( |A.1|) we write 



R(t) = ^t + V q t(1 - e' t/T ) + - f dt' f e-V-*">l T ttf)df. 

T JO JO 



(21) 



Therefore, 



E 



(R(t)-R(t)) 2 =2K 11 (t) 



2a ~ ft dt' f dt'e-V-*'^ fdt x W-e-^'A EUt^t") 



T JO 



where Kn(t) is given by Eq. @. Taking into account that 

E[amt)]=s(t'-t), 

we have 



E 



(R(t)-R(t)) 2 =2K u (t) 



2a 2 ft 



T JO 



dt' / dfe-V-*'^ \l - e-^"^ 



However, (see Eq. 



- Tdt' I" dfe-V-*"^ h - e -(«")/rl = Ki1 m 

T JO JO 1 J 



2 r t 



Hence, 



E 



(R(t)-R(t)f 



0. 



and R(t) is equal to R(t) in mean square sense. 
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3.3. The projected process when the initial velocity is in the stationary regime 

As we have mentioned, we are mainly interested in representing the asset dynamics when 
the initial velocity Vq is random and distributed according to the stationary pdf fllPP . 
We have shown in Section |2| that this basically implies the replacement of Ku(t) by 
K,(t). In such a case, the SDE for R(t) readsfjj] 



dR(t) = fidt + ^Jk(T - t)dW(t), 

where n(t) is given by Eq. ( |TTD and the dot denotes time derivative, that is 

k(t) =a 2 (l- e-*/ T ) . (22) 

We need the Ito lemma given in Appendix B for deriving the SDE for the stock S. Thus, 
according to Eq. (p.6| ), the effective dynamics for S = S e R is 

dS(t) 



S{f) \» + k{T - t)/2]dt + ^k{T - t)dW(t). (23) 

In this way, we have projected the two-dimensional O-U process (S, V) onto a one- 
dimensional price process which is a Wiener process with time varying drift and 
volatility. We also note that we need to specify the final condition of the process because 
the volatility \fk is a function of the time to maturity T — t, and this implies that the 
projected asset model depends on each particular contract. 

4. The option price on the projected process 

In this section we will present a generalization of the Black-Scholes theory assuming 
that underlying price is driven by the O-U process. We therefore eliminate the efficient 
market hypothesis but retain the other two requirements of the original B-S theory: the 
absence of arbitrage and the existence of a riskless strategy. 

We invoke the standard theoretical restrictions -continuos trading without 
transaction costs and dividends- and apply the original B-S method taking into account 
that the underlying asset is not driven by white noise but by colored noise modelled as 
an O-U process. 

The starting point of B-S option pricing is a portfolio which contains certain 
amounts of shares, calls and bonds. In this context, B-S hedging is only able to remove 
the call risk that comes from stock fluctuations. Therefore, we need to start from the 
effective one-dimensional market dynamics given by Eq. fl2~3"[) since otherwise we would 
not be able to remove risk fluctuations arising from dW(t). These fluctuations are 



only explicitly given in the projected SDE for the stock (see Section |0| for a deeper 
discussion on this point). 



ff Since R(t) and R(t) are equal in mean square sense we will drop the bar on R as long as there is no 
confusion. Thus, we will use R for the projected process as well. 
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Black-Scholes option pricing with the equivalent one- dimensional SDE 

As we have proved in Section |3], there exists an effective one-dimensional diffusion which 
describes the O-U process (@)-(|3|). Assuming that the effective one-dimensional price 
dynamics is given by Eq. (p3|), it is quite straightforward to derive the European call 
option price within the original B-S method. Following Merton (1973b) we define a 
portfolio compounded by a certain amount A of shares at price S, a quantity of bonds 
$, and a number \1/ of calls with price C, maturity time T and strike price K. We 
assume that short-selling is allowed and thus the value P of the portfolio is written 

P = - AS - $5, (24) 

where the bond price B evolves according to the risk-free interest rate ratio r. That is 

dB = rBdt. (25) 

The portfolio is required to obey the net-zero investment hypothesis, which means P = 
for any time t (Merton (1973b)). Hence, 

C = 5S + <j)B, (26) 

where 5 = A/ty and = $/\& are, respectively, the number of shares per call and the 
number of bonds per call. Due to the nonanticipating character of 5 and <ft we have 
(Bj6rk (1998)) 

dC = 5dS + <j>dB. (27) 

On the other hand, assuming that the market dynamics is described by Eq. (|23f) , the 
differential of the call also reads 

dC(S, t) = C t dt + C s dS + ^k{T - t)S 2 C ss dt, (28) 



where we have used the Ito lemma as expressed by Eq. flB.7j) of the Appendix B. From 
Eqs. (|3)-(|D and (||) we get 

C t + ^k{T - t)S 2 C ss + r5S -rC dt = [5 - C s }dS. 

Now the B-S delta hedging, 5 = C$, removes any random uncertainty in the option 
price. The partial differential equation for C(S,t) then reads 

C t = rC- rSC s - ~/i(T - t)S 2 C ss . (29) 

We note that the delta hedging is able to remove risk because we have projected the 
two-dimensional SDE (§)-(§) onto the one-dimensional process. In this way, we directly 
relate the differential of the stock dS(t) to the random fluctuations of the Wiener process 
dW(t) (see Eq. (|23|)). Without this projection, the B-S hedging is useless and the 
random fluctuations persist in the B-S portfolio. We will further discuss this situation 



in Section 5.1 
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4-2. The price of the European call 

For the European call, Eq. fl29|) has to be solved with the following "final condition" at 
maturity time T 

C(3,T) = max[S(T)-K,0], (30) 

where S(T) is the underlying price at maturity and K is the strike price. The solution 
to Eq. (p9|) subject to Eq. ( 130|) is a type of solution perfectly known in the literature 
(see, for instance, Hull (2000)). Thus, our final price is 

C ou (S, t) = S N(d° u ) - Ke~ r{T - l) N(d° u ), (31) 

where 

N(z) = (l/V^) f e~ x2/2 dx 



oo 



is the probability integral, and 



j OU _ln(S/K)+r(T-t) + K (T-t)/2 



d? u = v 1 ' \ ' v ^, (32) 

\I<T - t) 



<% U = d° u -y/ K (T-t), (33) 

with K,(t) given by Eq. (pi]). 

Equation ([31]) constitutes the key result of the paper. Note that, when r = 0, the 
variance becomes n(t) = a 2 t and the price in Eq. (|3~ID reduces to the Black-Scholes 
price: 

C BS (S, t)=S N(df s ) - Ke-< T -^ N(d^ s ), (34) 

where (iff have the form of Eqs. (|32|)-(|33D with K(T—t) replaced by cr 2 (T— t). Therefore, 
the O-U price in Eq. (^l|) has the same functional form as B-S price in Eq. when 
a 2 t is replaced by n(t). 

In the opposite case, r = oo, where there is no random noise but a deterministic 
and constant driving force (in our case it is zero), Eq. (|3l| ) reduces to the deterministic 
price 

C d (S, t) = max [S - Ke- riT - t] , o] . (35) 

We will now prove that Cou is an intermediate price between B-S price and the 
deterministic price (see Fig. 0) 

C d {S,t) < C OU (S,t) < C BS (S,t), (36) 

for all S and < t < T. In order to prove this it suffices to show that Cou is a monotone 
decreasing function of the correlation time r, since in such 

C OU (r = oo) < C ou {r) < C OU (r = 0). 

However, Couij — °°) = an d Couij = 0) = Cbs, which leads to Eq. (|3"B]). Let us 
thus show that Cou is a decreasing function of r for < t < T and all S. Define a 
function a as the derivative 

dCou , „\ 

a = ~TT- (37) 
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Figure 2. Relative call price C/S as a function of S/K for a given time to expiration 
T — t = 5 days. The solid line represents the O-U call price with r = 1 day and the 
dashed line is the B-S price. The dotted line is the deterministic price. In this figure 
the annual risk-free interest rate r = 5%, and the annual volatility a = 30%. 



Since the r dependence in Cou is a consequence of the variance n(t, r), we have 

a d K (T-t,r) „ 

a ~ 2«(T-t,r) ¥r VoC/ ' 

where Vou — dCou/da (see Section 0). But 

9<T dr t; T) = -a 2 [l - (1 + (T - t)/r)e-( T -^] < 0, 

for < t < T which is seen to be non positive. From Eq. fl67D below we see that Vou ^ 
for all S and < t < T. Hence, a < which proves Eq. ©. In Fig. | we plot the 
option price C as a function of the correlation time r and for three different values of 
the moneyness S/K. This figure clearly shows that C is a monotone decreasing function 
of r. 

Therefore, the assumption of uncorrected underlying assets (B-S case) overprices 
any call option. This confirms the intuition understanding that correlation implies more 
predictability and therefore less risk and, finally, a lower price for the option. In fact, 
we can easily quantify this overprice by evaluating the relative difference 

D = (Cbs — Cou)/Cbs- 

Figure § shows the ratio D(S,t), for a fixed time to expiration, plotted as a function 
of the moneyness, S/K, and for different values of correlation time r. We see there 
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Figure 3. Relative call price C/S as a function of r for a given time to expiration 
T—t = 10 days. The solid line represents the call price with S/K = 1 (ATM case). The 
dotted line is the call price when S/K — 1.03 (ITM case). The dashed line represents 
an OTM case when S/K = 0.97. We clearly see that C is a monotone decreasing 
function of t having its maximum value when r = (B-S case) and its minimum 
when t — > oo (deterministic price). The annual risk- free interest rate and the annual 
volatility are as in Fig. 2. 



that the ratio D is very sensitive to whether the call is in the money (ITM), out of the 
money (OTM) or at the money (ATM). The biggest difference between prices occurs in 
the case of OTM options. This is true because when S/K < 1, both Cbs an d Cou are 
small but Cbs Cou (see Fig. |). Depending on the value of correlation time r this 
implies that D is approximately equal to 1. 

Another interesting point is the behavior of D as a function of the expiration time 
T — t. In this case, D behaves quite differently depending on whether the call is in, 
out, or at the money. This behavior is evident in Figs. |5| and [| Figure [5] shows D(S, t) 
as a function of expiration time T — t for an OTM option (S/K = 0.95) and the ATM 
option (S/K = 1.00) and for two different values (1 and 5 days) of the correlation time. 
Note that B-S notably overprices the option, particularly in the OTM case. In Fig. ^ 
we show plots of D(S,t) as a function of t for an ITM option (S/K = 1.05). This 
exhibits completely different behavior since the B-S overprice is considerably less (no 
more than 7%). Moreover, contrary to the ATM and OTM cases, the relative difference 
D(S, t) is a non monotone function of T — t, having a maximum value around one or two 
weeks before maturity. Although perhaps the most striking and interesting feature is the 
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Figure 4. D(S, t) is plotted as a function of S/K for T — t = 10 days and r = 1 day 
(solid line), r = 2 days (dashed line) and r = 5 days (dotted line). Other parameters 
used to generate the figure are r — 5% per annum and a — 30% per annum. 

persistence of the B-S overprice far from maturity regardless the value of the correlation 
time. This is clearly shown in Table [I] where we quantify the ratio D in percentages for 
different values of moneyness, time to expiration and correlation time. 

5. An alternative derivation of the call price 

In this section and the next, we present two different and alternative derivations of the 
final call price Cqu- The first of these derivations is based on an extension of the B-S 
theory but now starting from the two-dimensional diffusion (0)-([3P and with a different 
portfolio than the usual one. A second derivation, briefly outlined in the next section, 
uses the equivalent martingale measure method. Both derivations arrive at the price 
formula (0), thus showing the consistency of the pricing methods. 

We will first apply the original B-S method starting from the two-dimensional O-U 
process (!)-(§) instead of the equivalent process (|23D . Unfortunately, this procedure 
yields a trivial expression for the price of the option (see below) and is therefore useless. 
To avoid this difficulty we will define a different portfolio which is the first step towards 
the generalization of both B-S equation and formula. 
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Figure 5. D(S, t) is plotted as a function of T — t (in logarithmic scale) for fixed values 
of moneyness. The solid lines represent ATM options, the thick line corresponds to 
t = 1 day and the thin line corresponds to r — 5 days. The dashed lines represent an 
OTM option with S/K = 0.95, the thick line corresponds to r = 1 day and the thin 
line corresponds to t = 5 days (r and a as in Fig. 2). 



5.1. The Black-Scholes method for the two-dimensional O-U process 



We assume that market prices are driven by an O-U process as shown in Eqs. (@)-(0) 
and that the portfolio is given by Eq. fl24]) . That is, C = 5S + (pB and 

dC = 5dS + <j>dB. 

Let us now apply the original B-S method starting from the two-dimensional O-U 
process (@)-@ instead of the equivalent process (p3|). Using the Ito lemma for a singular 
two-dimensional diffusion (see Appendix B), 

r 2 



a 



dC(S, V, t) = C s dS + CydV + C t dt + —C vv dt, 

_ / 

and taking Eqs. (|25| ) and ( p7[ ) into account, we write 



(38) 



a 



C t + —C vv -r(C-S5) 
It 



dt + (C s - 5)dS + CydV = 0. 



Now the assumption of delta hedging 5 = Cs, turns this equation into 



a 



C t - r(C - SC S ) + —C vv 
It 



dt + C v dV = 0. 



(39) 
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Figure 6. D(S, t) is plotted as a function of the expiration time T — t (in logarithmic 
scale) for an ITM option with S/K = 1.05. The solid line corresponds to r = 1 day 
and the dashed line to r = 5 days (r and a as in Fig. 2). 



Equation (|39| ) is still random due to the term with dV representing velocity fluctuations 
(see Eq. (||)). In consequence, Black-Scholes delta hedging is incomplete since it is not 
able to remove risk. In this situation, the only way to derive a risk-free partial differential 
equation for the call price is to assume that the call is independent of velocity Then, 
C v = and Eq. © yields 

C t + rSC s -rC = Q. (40) 

According to the final condition for the European call, C(S,T) = max[5(T) — A', 0], 
the call price is C(S, t) = max S — Ke~ r<yT ~ l \ . Note that this is a useless expression 
because it gives a price for the option as if the underlying asset would have evolved 
deterministically like the risk-free bond without pricing the random evolution of the 
stock. In fact, there is no hint of randomness, measured by the volatility a, in Eq. fllPP . 

The main reason for the failure of B-S theory is the inappropriateness of B-S hedging 
for two-dimensional processes such as O-U price process (@)-(||)@. Indeed, delta hedging 
presumably diversifies away the risk associated with the differential of asset price dS(t) 
given by Eq. Nevertheless, what we have to hedge is the risk associated with dV(t) 
given by Eq. (|3|), which contains the only source of randomness: the differential of the 
Wiener process dW(t). All of this clearly shows the uselessness of the B-S delta hedging 
for the two-dimensional O-U process. Note that we must relate in a direct way the 



| A similar situation appears in the stochastic volatility models (Scott (1987)). 
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Table 1. Relative call price differences in percentages. Values of D x 100, where 
D = (Cbs — Gou)/Cbs- T — t is the expiration time in days. Correlation times r are 
1, 2, and 5 days. The rest of columns are divided in three blocks corresponding to a 
different values of the moneyness S/K . From left to right blocks represent the OTM, 
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differential dS(t) with the random differential dW(t), otherwise we will not be able to 



remove risk. This is indeed the case of the projected process ( |23"D which leads to the 
European call price, Eq. ([H]). However, if we do not want to project the process and 
maintain the two-dimensional formulation we have to evaluate the option price 

from a different portfolio. We will do it next by defining a new portfolio which will 
allow us to preserve the complet market hypothesis and remove the random component 
dW(t). 

5.2. The option pricing method with a modified portfolio 



We present a new portfolio in a complete but not efficient market. The market is 
still assumed to be complete, in other words, there exists a portfolio with assets to 
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eliminate financial risk. However, we relax the efficient market hypothesis by including 
the correlated O-U process as noise for the underlying price dynamics. 

Now, our portfolio is compounded by a number of calls \I/ with maturity T and 
strike K, a quantity of bonds $, and another number of "secondary calls" on the 
same asset, but with a different strike K' and, eventually, different payoff or maturity 
time. Note that in the new portfolio there are no shares of the underlying asset. Thus, 
instead of Eq. (|^f) , we have 

p = mc - m'C - <&b. (41) 

After assuming the net-zero investment, we obtain 

C = (j)B + ipC, (42) 

where = is the number of bonds per call, and ip = ^'/^ is the number of 

secondary calls per call. We proceed as before, thus the nonaticipating character of (j> 
and i[) allows us to write 

dC = <j>dB + i)dC (43) 

and, after using Ito lemma (581) for both dC and dC, some simple manipulations yield 



a 2 



C t + —C vv -rC + (fi + V)SC s 
It 



a 2 



^C' vv -rC ! + ^ + V)SC' & 



dt = tyC'y - C v ) dV. (44) 



This equation can be transformed to a deterministic one by equating to zero the term 
multiplying the random differential dV(t) given by Eq. (|3|). This, in turn, will determine 
the investor strategy giving the relative number of secondary calls to be held. Thus, 
instead of B-S delta hedging, we will have the "psi hedging" : 

V> = % (45) 



Then 



a 



V 



a 2 



1 C t + —C vv -rC + {^ + V)SC s 
It 



Cv 



o~ 2 



C t + Y Cyy - rC + (/X + V)SCg\ . (46) 



This equation proves, as otherwise expected, that the call has the same partial 
differential equation independent of its maturity and strike. This has been suggested in 
a more theoretical setting for any derivative on the same asset (Bjork (1998)). 

On the other hand, the two options C and C have different strikes. Then, 
analogously to the separation of variable method used in mathematics (Mynt-U (1987)) 
and proceeding in a similar way to that used in the study of SV cases, both sides of 
Eq. (|46D are assumed to be equal to an unknown function X(S, V, t) of the independent 
variables S, V, and t. We thus have 

9 

C t + °-C vv + (fi + V)SC S -rC = \C V . (47) 
It 



Option pricing and perfect hedging on correlated stocks 



19 



In the stochastic volatility literature, the arbitrary function X(S,V,t) is known as the 
"risk premium" associated, in our case, with the return velocity (Scott (1987); Heston 
(1993)). In the Appendix C we show that the risk premium A is given by 

X(S,V,t) = -. (48) 
r 



A substitution of Eq. fl48|) into Eq. ([47]) yields a closed partial differential for the call 
price C(S, V, t) which is 

C t + t-C yv - -C v + (// + V)SC S - rC = 0. (49) 
It t 

For the European call, Eq. ([49]) has to be solved with the "final condition" ([30]) at 
maturity time which is C(S, V, T) = max[5(T) — K, 0]. The solution to Eq. ([49]) subject 
to Eq. ([30]) is given in Appendix D and reads 

C(3, V, t) = e- r{T - t] [Se^ T - t ' v) N{z l ) - KN(z 2 )] , (50) 

where Z\ = Zi(S, V,T — t), z 2 = z 2 (S, V,T — t) are given by Eq. ( |L).10| ) of Appendix D, 
and 

(3(t,V) = m(t,V) + K n (t)/2, 

where m(t, V) and Knit) are given by Eqs. (§) and @. 

The option price ( |50"|) depends on both the price S and the velocity V of the 
underlying asset at time t, i.e., at the time at which the call is bought. They are 
therefore the initial variables of the problem. However, while the initial price S is 
always known, the initial velocity V is unknown. The velocity is thus assumed to be in 
the stationary regime so that its probability density function is as shown in Eq. (|10|). 
We therefore average over the unknown initial velocity and define C by 



C(S,t)= / C(S,V,t)p st (V)dV, (51) 

and from Eqs. ( |TUD and flSO]) we have 

C(S, t) = e- r(T -* } [Se^'^Nizt) - KN(z 2 )] , (52) 

where 

P(t) =id + n(t)/2, (53) 
n(t) is the variance defined by Eq. (pTT]) , and Z\ t2 are given by Eq. ( p.ll| ) of Appendix D. 



As mentioned above, Eq. (15*2] ) cannot be our final price yet because it still depends on 
the mean return rate \i. This rate could differ depending on whether \i is estimated by 
the seller or buyer of the option and thereofore, in Eq. ([52]), there are hidden arbitrage 
opportunities. 

Therefore, we must proceed in a similar way as in the martingale option pricing 
theory of Eq. ([62]) and define the final call price, Cou(S,t), as price C when (3{t) is 
replaced by rt. That is: 

C OU (S,t)=~C(S,t)\ m ^ t , (54) 
and this price completely agrees with the one derived in Section |3| (see Eq. ([31]) ). 
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Suppose we start from the modified portfolio ( pE2|) but assuming that the share price is 
given by the projected process (|^) instead of the two-dimensional O-U process 
In this case, one can obtain the same option price as before (cf. Eq. (|31~D ). However, 
the hedging strategy will be given by the following function 

C s 



C' 



(55) 



Let us prove this. We start from Eq. (fi3|): 

dC = <j)dB + tpdC, 
Now, instead of Eq. we have (see Ito lemma (|28] 

(C t + h(T-t)C ss -rC 



i,{C , t + -k(T-t)C 88 -rC 
= tyC' s -C s )dS. 



dt 



(56) 

And the removal of risk implies Eq. (|55p . The psi hedging given by Eq. (|53|) is equivalent 
to the psi hedging defined in Eq. (fi5| ) although now it is represented in terms of the 
final price Cou(S,t) instead of the intermediate price C(S, V,t). Substituting Eq. (p5|) 
into Eq. (^) and reasoning along the same lines as above (see Eq. ($7\)) we obtain 



C t + -k(T-t)S 2 C ss -rC 



AC, 



(57) 



where A = X(S, t) is the "risk premium" for the effective process which is now obviously 
independent of the velocity V. Combining Eqs. (p3|), (p8|) and (|57|), we get 



(5f 



dC(S, t) = [rC + [- + n + ]-k{T - t)] SC s )dt 



+ yJk(T-t)SC s dW(t). 
Hence, the conditional expected value of dC reads 



E[dC\C] 



rC 



A 1 /m X 

- + M+2«(T-t) 



SC S dt, 



but the equilibrium of the market implies that i?[dC|C] = rCdt. Therefore, 



(59) 



(60) 



A = -S 
and Eq. fl5"7|) reads 

C t + l -k(T - t)S 2 C ss - rC + [/i + l -k(T - t)] SC S = 0, 
Finally, the absence of arbitrage opportunities requires the replacement (see Eq. fl53|)) 

H + k(T-t)/2 — ► r. 
Thus, the option price equation is 



C t = rC- rSC s - \k{T - t)S 2 C ss 



(61) 
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which agrees with Eq. (|29|). 

Note that both procedures, the original B-S method presented in Section |3| and our 
method, result in the same partial differential equation for the call price. However, each 
method uses a different hedging strategy because they start from a different portfolio. 



6. The call price by the equivalent martingale measure method 

As was shown by Harrison and Kreps (1979) and Harrison and Pliska (1981), the B-S 
option price can also be found using martingale methods. This is a shorter, although 
more abstract way, to derive an expression for the call price. The main advantage is 
that one only needs to know the probability density function governing market evolution 
which, in turn, allows one to obtain an option price in situations where B-S assumptions 
are not applicable. The drawback is that one is not sure of whether the price obtained 
by martingale methods is the fair price of the call because of the omission of arbitrage 
and hedging. 

We will now show that, in the present case, the price obtained by martingale 



methods completely agrees with our extended B-S price (0). The equivalent martingale 
measure theory imposes the condition that, in a "risk-neutral world", the stock price 
S(t) evolves, on average, as a riskless bond (Harrison and Pliska (1981)). 

Let p*(S,t\S ,t ) be the equivalent martingale measure associated with asset price 
S(t) conditioned on S(t ) = So- Define the martingale conditional expected value 

E*[S(t)\S ] = / Sp*(S,t\S ,t )dS. 

Jo 

Then the risk-neutral assumption requires that 

E* [S(t)\S }=S e r «- t °\ 
where r is the constant spot interest rate. On the other hand, 

E[S(t)\S ]= / Sp{S,t\S ,t )dS. 
Jo 

Assuming that the initial velocity is in the statioanry regime, the marginal density 
p(S, t\St ) is given by Eq. ( |A.16|) of Appendix A. Therefore, 

E[3(t)\3 ] = S Q exp W-h)], 

where (3{t) = fit + «(t)/2 with n(t) given by Eq. ([□]). We thus see that the equivalent 
martingale measure is accomplished by the replacement 

P(t) — ► rt. (62) 

In consequence, 

i 

p*(S,t\S ,t ) 



£y27TK(t - t ) 

1 2 K (f-i„) /' * ' 
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which is the so called "the risk-neutral pdf ' for the stock price and it is a consequence 
of the absence of arbitrage demmand. Now, it is possible to express the price for the 
European call option by defining its value as the discounted expected gain due to holding 
the call. That is (Harrison and Pliska (1981)), 

C*{S,t) = e- r(T ^E*{max[S{T) - K,0]\S(t) = S} 

POO 

= e - r(T -* ) / (S' - K)p*(S', T\S, t)dS', (64) 
Jk 

and the final result for the call is obtained by calculating the expected value with the 
equivalent martingale measure defined in Eq. fl5"B|). The martingale price agrees exactly 
with our previous price in Eq. (|HT|), C*(S,t) = Cou(S,t). We can thus say that, in the 
O-U case, both option pricing methods are completely equivalent although martingale 
theory does not require the construction of a portfolio and ignores any hedging strategy. 



7. Greeks and Hedging 



We briefly derive the Greeks for the O-U case. Since the O-U call price has the same 
functional form as the B-S price but replaces er 2 (T — t) by k(T — t), the O-U Greeks 
will have the same functional form as B-S Greeks with the same replacement except 
for Vega, V = dC/da, and 9 = dC/dt. Thus, for 5 = OC/dS, 7 = d 2 C/dS 2 , and 
p = dC/dr, we have (Hull (2000)) 

p -(d? u ) 2 /2 

N(d° u ) 



ou 



lou 



S^2<kk{T - 1) 

pou = K(T - t)e~ r ^N(d° u ). (65) 

Since d®^ > dfi f° r an $ an d t and N(z) is a monotone increasing function, we see 
that Sou > 8bs and pou > Pbs- Hence, the O-U call price is more sensitive to changes 
in stock price and interest rate than the B-S price. 

On the other hand, from Eq. ([H]) and taking into account the identity 

(66) 



we have 



and 



SW'(di) - Ke~ r{ - T ~ t) N\d 2 ) = 0, 

(S/aMT-t)^} 1 ^-^ 2 / 2 , 



ou 



0, 



ou 



_tf e -»-(T-t) 



rN(d 



a' 



l_ e -{T-t)M 



-(dOUy/2 



(67) 



(68) 



2^2ttk(T - 1) 

bmce d° u > df s , one can easily see that Vou < Vbs f° r an values of S/K, T — t and 
t. Thus our correlated call price is less sensitive to any change of underlying volatility 
a than is the B-S price. 

We conclude with the psi hedging. For the two-dimensional O-U case the hedging 
strategy is given by the function i/;(S,V,t) specifying the number of secondary calls 
to be hold. However, the hedging given by Eq. (EE]) depends on the velocity V and 
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Figure 7. Hedging in terms of the moneyness. Psi hedging and delta hedging as a 
function of the moneyness S/K. The solid line represents psi hedging when r = 1 day, 
the time to expiration is T — t = 20 days, and the exercising price of the secondary 
call is K 1 — 0.9K. The dotted line corresponds to the delta hedging still assuming 
the O-U asset model with the same correlation and expiration time. The dashed line 
corresponds to B-S delta hedging (r and a as in Fig. 2). 



is not expressed in terms of the final call price Cou — C(S,t). As we have shown in 
Section |5.3| , psi hedging in terms of Cou can only be derived from the effective one- 
dimensional process fl2"3"|). In this case, the removal of the randomness coming from dS 
implies that hedging is given by Eq. (|55|). Since Cs = Sou, we see from Eqs. (|55|) and 
Eq. (HP that 

Now, we take the secondary option to be an European call with maturity T and 
exercising price K' < K, where T and K refer to the primary option. We plot in Fig. [7| 
the psi hedging as a function of the moneyness. We see there that the ip hedging is 
always greater than Sou an d Sbs hedgings. Since N(d[) — > 1 when K' — » 0, the psi 
hedging approaches to the delta hedging Sou as the moneyness of the secondary call 
tends to infinity. This is consistent with the fact that secondary calls have the same 
price as the underlying stock when its exercising price is zero , i.e. C —>■ S as K' — > 
(see Eq. (|3~l"D). Therefore, having secondary calls with exercising price equal to zero is 
equivalent to own underlying shares and the O-U psi hedging coincides with the O-U 
delta hedging. 
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As we have mentioned, psi hedging ip indicates the number of secondary calls per call 
to be hold if we follow a risk-free strategy with the modified portfolio (fjip. Therefore, 
the money invested to carry out this strategy is given by tpC . That is 

1>C = |[|| [SN(d[) - K'e-^N(d' 2 )} , (70) 

where we have combined the Eqs. ([$!]) and (|69|) . On the other hand, delta hedging also 
indicates the number of shares per call to be hold in a risk- free strategy with the B-S 
portfolio (p4[) . And, analogously, the money necessary to perform this strategy is 

5S = SN(di), (71) 

where 5 is given by Eq. (^). We compare these quantities in order to know which 
hedging is cheaper for the investor. From Eqs. (|70|) and ([H]), we see 

i)C :j K'e- r ^ N(d' 2 ) 



butg 



SS S N(d[) 



K ' e -r(T-t) N , d n 

< — < 1. 

S N{d[) 



Therefore, ipC < 5S and psi hedging is always less expensive than delta hedging. 
Note that when K' — > both strategies have the same cost. In Fig. ||] we plot, as a 
function of moneyness, the relative psi hedging cost, ipC'/K, along with the relative 
delta hedging cost, SS/K. We see there that psi hedging is considerably less expensive 
than delta hedging and this difference increases with moneyness. Indeed, for an ATM 
call (S/K = 1.00) and with parameter values as that of Fig. ||, delta hedging is 
approximately 800% more expensive than psi hedging. 

Combining Eqs. (0), (^) and (|66|) one can easily show that when r = the O-U 
psi hedging is i) BS = N(df s )/N(df s ') I, where the prime refers to the secondary call. 
Since 5 = C$ = N(di), we have 



BS 



Finally, for the secondary call, whose exercising price goes to zero, 5' BS — > 1 and, again, 
B-S psi hedging and B-S delta hedging coincide. 



8. Conclusions 



We have developed option pricing with perfect hedging in an inefficient market model. 
The inefficiency of the market is related to the fact that the underlying price variations 
are autocorrelated over an arbitrary time period r. In order to take these correlations 
into account we have modelled the underlying price S(t) as a singular diffusion process 

§ This is straightforward to prove from Eq. since ipC > 0. 

| We use the subscript BS in ipgs to indicate that this hedging refers to an uncorrelated stock (r = 0), 
as in the B-S world 
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moneyness 



Figure 8. Relative hedging costs ipC /K and SS/K as a function of the moneyness 
S/K. The solid line represent psi hedging cost when r = 1 day, the time to expiration 
is T — t = 20 days, and the exercising price of the secondary call is K' — Q.9K. The 
dotted line corresponds to the delta hedging with r = 1 day and T — t = 20 days (r 
and a as in Fig. 2). 

in two dimensions (O-U process) instead of the standard assumption that S(t) is a one- 
dimensional diffusion given by the geometric Brownian motion with constant volatility. 

The option pricing method has been developed by keeping perfect hedging with a 
riskless strategy which finally results in a closed and exact expression for the European 
call. Our pricing formula has the same functional form as the B-S price but replaces the 
variance of the Wiener process by the variance of the O-U process. The O-U variance, 
n(t), is smaller than the B-S variance, a 2 t, which implies that the equivalent volatility 
in the O-U case is lower than B-S volatility^]. But less volatility implies a lower option 
price. We have indeed proved that the B-S call price is always greater than the O-U 
price. In other words, the assumption of uncorrelated assets overprices the European 
call. This agrees with the fact that correlation, which can be regarded as a form of 
predictability, implies less risk and therefore a lower price for the option. We have 
quantified this overprice and showed that B-S formula notably overprices options and, 
more strikingly, that the overprice persists for a long time regardless of the strength of 
correlations. We have also analyzed the sensitivity of the O-U price to several conditions. 

% Since the volatility a is the square root of the variance per unit time, one can define, in the O-U 
case, an equivalent volatility by aou = \J where the dot denotes time derivative. From Eq. ( |TT| ) 
we see that ctqu /& = vl — e~*/ r < 1. 
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Thus we have proved that while Cou is more sensitive to changes in the interest rate 
and stock price than Cbs, h is also less sensitive to any change of the volatility. The 
practical consequences of this are nontrivial. 

The option price and the hedging strategy have been obtained using two different 
approaches. The most straightforward way of getting the call price is by means of a 
projection onto a one-dimensional process with a time-varying volatility. A second way 
of obtaining the option price starts with the complete two-dimensional O-U process (fj)- 
(§). This is a longer procedure but opens the door to a new hedging strategy: the psi 
hedging. We have therefore two ways of acheiving the perfect hedging: the usual one 
consisting in holding underlying assets (delta hedging), and the second one which uses 
secondary callsQ instead of assets (psi hedging) . We have shown that this last strategy 
can be considerable less expensive than the delta hedging and can avoid a possible lack 
of liquidity of underlying shares. Finally, the proportion of secondary calls to be held, 
i.e., the psi hedging, converges towards O-U delta hedging when the exercising price of 
the secondary call tends to zero. 

In practice our method of valuation requires the estimate of one more parameter, 
the correlation time, than in the B-S Wiener case. Assuming that the underlying asset 
is driven by O-U noise one can find an estimate for the correlation time r by evaluating 
the variance n(t) of the asset return. Once one has an estimate of this variance the 
correlation time is given in Eq. flTTl). 

We finally mention that one interesting extension of the valuation method presented 
is to the American option. Although this case is more involved, one is probably able 
to obtain, at least an approximate or a numerical result using a combination of first 
passage times and martingale methods, as recently presented by Bunch and Johnson 
(2000). In any case we believe that the effects of autocorrelations on the valuation of 
an American option will be even more critical than for the European call. This case is 
under present investigation. 
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+ The construction of the portfolio with secondary calls is one simple way of proceeding. Obviously, 
any other secondary derivative on the same asset would serve. 
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Appendix A. Mathematical properties of the model 

We present some of the most important properties of the model given by the pair of 
stochastic equations in Eqs. (||) and (|J). Their formal solutions are 

V(t) = V e- (t - to)/T + - /V(*-*'}/ W(i'), 

and 

R(t) = fi(t - t ) + V t(1 - e -(*-*°>/<-) 

+ <L f dt ' f e-V-^dWit"), (A.l) 

T Jto Jto 

where we have assumed that the process begun at time to with initial velocity Vq and 
return R = 0. The return R(t) has the following conditional mean value 

E[R(t)\V ] = fi(t - to) + r(l - e-^^Vo, 

and variance 

Var[R(t)|V ] = o- 2 [(t - t ) - 2r (l - e'^ ^) + T - (l - e" 2 ^-^)" . 

Since {R{t), V(t)) is a diffusion process in two dimensions, its joint density p(R, V, t) 
satisfies the following Fokker-Planck equation (Gardiner (1985)) 

V o 2 

p t = -(/i + V)p R + — p v + j-^Pvv- (A.2) 

This is to be solved subject to the initial conditions R(t ) = and V(t ) = V , that is 

p(R,V,t \V ,t ) = 5(R)5(V -V ). (A.3) 

A first step towards solving the problem ( |A.2[ )-( |A^D is the definition of the joint 
Fourier transform 

/oo /*oo 
dRe iaR dVe i(3v p{R,V,t). 
-oo J —oo 

Then problem ( |A.2| )-( |A7^ ) becomes 

d t p = iafip +{a- f3/r)adpp - {a 2 /2r 2 )(3 2 p, (A.4) 

p{a,P,t = 0) = e im . (A.5) 

We look for a solution of the form 

p(a,/3,t) = exp{i[ami(t) + /3m 2 (t)]} 

x exp{-[^ n (t)a 2 + K 12 (t)af3 + K 22 (t)(3 2 ]/2}, (A.6) 

where mj(t) and Kij{t) are functions to be determined. We substitute Eq. (|A.6| ) 



into (|A.4j) and identify term by term. We have 

mi = fi + m 2 , m 2 = -m 2 /T; 

K 22 + (2/t)K 22 = a 2 r 2 , K l2 + (l/r)K 12 = 2K 22 (t), K u = 2K U , 
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with the intial conditions, according to Eqs. ( |A.5| )-( [7C6D , given by 

m 2 (0) = V (h m 1 (0) = K ij (0) = = 1,2). 

The solution reads 
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mi(t) = fit + V t (l - e- t/T ) 
and Kij{t) are given by 



m 2 



(t) = V e-^, 



K n {t) = a 2 



t - 2r (l 



-t/r 



1 - e 



-It It 



K 



12 



(t) = £ (1 - e-"i 



K 22 (t) 



a" 
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1-e 



-2t/r 



(A.7) 
(A.8) 



The inverse Fourier transform of Eq. 
p{R,V,t\V ,t ) 



yields the Gaussian density 



2irJdet[K(t - *„)] 



cxp < 



2if 22 (t - 1 ) 



^ 22 (t - to)(i? - m{t - t , Vo)) - K n (t - t ) (V - V^-^) 



2K 22 (t-t )det[K(t-t )} 



where 



and 



(A.9) 
(A.10) 
(A.ll) 

Notice that the joint density ( |A.9| ) is a function of the time differences t — to where to is 
the initial observation time, so that the two-dimensional diffusion (S(t), V(t)) is a time 
homogeneous process and, without loss of generality, we may assume that to = 0. 
The marginal pdf of the velocity V(t), 



det[K(t)] = K n {t)K 22 {t) - K 2 12 {t). 



m(t, Vq) = /it + V t (l - e- t/T ) . 



/oo 
p(R 7 V,t\V )dR, 
-oo 



IS 



p{V,t\V ] 



2nK 22 (t) 



exp 



(V-Voe-^y 
2K 22 (t) 



(A.12) 



In the stationary regime (t — > oo) we find a normal density independent of the initial 
velocity: 

1 



Pst(V) 



-TV 2 /a 2 



'7r(cr 2 /r) 

Analogously, the marginal density of the return R(t), 

/oo 
p(R,V,t\V )dV, 
-oo 



(A.13) 
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is 



p(R,t\V ) 



exp 



[R-m(t,V )f 
2K n (t) 



(A.14) 



/2irK n (t) 

If we assume that the initial velocity Vq = V(0) is a random variable distributed 
according to the pdf in Eq. ( A.13Q . We can therefore average the above densities to 
obtain a pdf independent of Vq. That is, 



p(R,V,t) 



p(R,V,t\V ) Pst (V )dV , 



and similarly for the marginal pdf's p(R,t) and p(V,t). Since we are mainly interested 
on the marginal distribution of the return we will give its explicit expression. Thus, 
from Eqs. (|A.13 ) and ( |A.14 ) we have 



p(R,t) 



27TK(t) 



: exp 



{R - fdf 

2K(t) 



(A.15) 



where K(t) is given by Eq. (pTTj) . Alternatively, the distribution of the underlying price 
S = Soe R is given by the log-normal density 



p(S,t\S 



SJ2lTK(t) 



exp 



{lnS/S -fity 
2n(t) 



(A.16) 



From this we easily see that the conditional probability p(S',T\S,t) when t < T is 



p(S',T\S,t) 



S>J2hk(T - t) 



: exp 



[In S'/S - /i(T ~t)f 
2k(T - 1) 



(A.17) 



Appendix B. The Ito formula for processes driven by O-U noise 

In this Appendix we generalize the Ito formula for processes driven by Ornstein- 
Uhlenbeck noise. This is applied to the share price S(t) which is governed by the 
pair of stochastic equations (0)-(@) 



dS(t) = S(fi + V)dt, dV(t) = -—dt + -dW. 

T T 



(B.l) 



Consider a generic function f(S,V,t) which depends on all of the variables that 
characterize the underlying asset. The differential of f(S, V,t) is defined by 

df(S, V, t) ee f(S(t + dt), V(t + dt),t + dt) - f(S(t), V(t),t). (B.2) 

But the Taylor expansion of (|B.2|) yields 



(B.3) 



df(S,V,t) = f s dS + f v dV + f t dt 



+ \fssdS 2 + \f vv dV 2 + fsvdSdV + 



where the expansion also involves higher order differentials such as (dt) 2 , (dS) 3 , (dV) 3 , 
etc. However, the differential of the Wiener process, dW, satisfies the well-known 
property, in the mean-square sense, dW(t) 2 = dt (Gardiner (1985)). And from the 
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pair of equations ( |B.1| ) we then see that dS 2 is of order dt 2 while dV 2 is of order dt and 
dSdV is of order dt 3 / 2 . Therefore, up to order dt, Eq. ( |B.3| ) reads 

■2 



df(S, V, t) = f s dS + f v dV + f t dt + ^fvvdt, 



(B.4) 



which is the Ito formula for our singular two-dimensional process ( 
Suppose now we start from the effective one-dimensional SDE (|l 



dR{t) = ^dt + ^k{T - t)dW{t). (B.5) 

We will prove that the corresponding SDE for the stock price defined as S = Soe R is 
given by Eq. (p^). In effect, substituting Eq. ( |B.5| ) in the Taylor expansion 

dS(R) = S R dR + ^S RR dR 2 + ■ ■ ■ , 

neglecting orders higher than dt and taking into account that dR 2 = k(T — t)dt (in mean 
square sense), we finally obtain 



*W) = + k(T - t)/2) dt + y/k(T-t)dW(t), 



(B.6) 



which is Eq. (||). 

Moreover, we can also give the differential of a generic function f(S,t) when 
underlying obeys SDE ( p.6|) . In this case, we have 

df(S, t) = f s dS + f t dt + ~k(T - t)S 2 f ss dt, (B.7) 
where again we have neglected higher order contributions than dt. 



Appendix C. A derivation of the risk premium 

We proceed to find a closed expression for the arbitrary function A(5*, V, t) that appears 
in Eq. (0). The call price C is a function of S, V , and t. We now consider this function 
taking into account that S = S(t) and V = V(t) follow Eqs. (0) and (|), respectively. 
This therefore allows us to evaluate the random differential dC using the Ito lemma, as 
a result we find that 

r 2 



dC 



a 



C t + ( f i + V)SC s + —C vv 
Zt 



dt + C v dV. 



After using Eqs. (f|7]) and (|3|), we have 



dC 



dt + -C v dW. 

T 



The expected value of dC, on the assumption that C(t) = C is known, reads 



E[dC\C) 



rC + (A - -)C V 



dt. 



We claim that this average must grow at the same rate as the risk-free bond: 
E[dC\C\ = rCdt, 



(C.l) 

(C.2) 
(C.3) 
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since otherwise the option would not be in equilibrium (Hull (2000)). In some sense, this 
assumption is similar to that of the equivalent martingale measure demand expecting 
that markets grow in average as the risk- free bond (Harrison and Pliska (1981)). 

The substitution of Eq. ( |C3| ) into Eq. ( |U.2| ) yields the following expression for the 
risk premium X(S, V, t): 
V 



A 



T 



(C.4) 



Appendix D. Solution to the problem in Eqs. (|49|)-(|3~0|) 



We will solve Eq. (|49| ) subject to the final condition in Eq. 
independent variable Z 

S = e z , 



where the domain of Z is unrestricted. The problem posed in Eqs. 



rC-(fi+ V)C Z + -C v - °—C vv , 
t 2r 



Define a new 
(D.l) 

9|)-(|30D now reads 
(D.2) 
(D.3) 



C(Z,V,T) = max[e z -K, 0]. 
The solution to this problem can be written in the form 

/oo /*oo 
dZ' dV'max[e z - K]G(Z,V,t\Z',V',T), (D.4) 
-oo J —oo 

where G(Z,V,t\Z',V',T) is the Green function for the problem (Mynt-U (1987)), i.e., 
G(Z, V, t\Z', V, T) is the solution to 



G t = rG-(fi + V)G Z + -G v 

T 



o- 

2t" 



G 



VVi 



with the final condition 

G(Z, V, T\Z', V, T) = 5(Z - Z')6(V - V) 



(D.5) 



(D.6) 



where 5(X — X') is the Dirac delta function. Define G = e rt G, then the final-value 
problem in Eqs. ( p.5| ) and ( p.6| ) reads 

_ _ V- a 2 - 

Gf — — (/i + V)Gz H Gy — ^— Gy\ . 



T 2t 

G{Z, V, T\Z', V, T) = e~ rT 5{Z - Z')5{V - V). 



(D.7) 
(D.8) 



Note that Eq. ( p.7| ) is the backward equation corresponding to Eq. ( |A.2|) . Therefore, 
Eq.( [A.9| ) permits us to write the solution to the problem posed in Eqs. (p.7| )- (|D.8j ) 
(Gardiner (1985)). This solution implies that G is 

G(Z, V, t\Z', V, T) = 1 expi -r(T - t) 



V' - Ve- {T - l ^ T 



27rJdet[K(T -t)} 



2K 22 (T-t) 



K 22 (T -t)(Z'-Z + m(T - t, V)) - K n {T - t) (V 1 - Ve 



(T-t)/r 



2K 22 (T-t) det[K(T-t)] 



(D.9) 
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where det[K(t)}, K^t), and m(t,V) are denned in Eqs. ( prT0|) -( pOTl) . 

Substituting Eq. ( p.9| ) into Eq. ( |JJ.4j ) and finally reverting to the original variables 
we obtain Eq. (pj) with 

\n(S/K) + m(T — t,V) + K n (T - t) 



zi = , z 2 = z x - JK U (T - t). 

/K u (T-t) 



(D.10) 

Finally it can be shown, after some lengthy but simple manipulations, that the functions 
zi,2 = zi,2(S,T — t) appearing in the averaged price C(S,t), Eq. (|52|) , are given by 

]n(S/K) + u(T — t) + k(T — t) i— r , 

zi = W ' 7 ; - L , z 2 = z 1 -J K (T-t), (D.ll) 

y/<T - 1) 

where n(t) is given in Eq. (|Tl|). 
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